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Abstract 

It was pointed out in a previous paper that although neither the Klein-Gordon equation nor the Dirac 
Hamiltonian produces sound solitary free-particle relativistic quantum mechanics, the natural square-root 
relativistic Hamiltonian for a nonzero-mass free particle does achieve this. Failures of the Klein-Gordon 
and Dirac theories are reviewed: the solitary Dirac free particle has, inter alia, an invariant speed well in 
excess of c and staggering spontaneous Compton acceleration, but no pathologies whatsoever arise from the 
square-root relativistic Hamiltonian. Dirac's key misapprehension of the underlying four-vector character 
of the time-dependent, configuration-representation Schrodinger equation for a solitary particle is laid bare, 
as is the invalidity of the standard "proof" that the nonrelativistic limit of the Dirac equation is the Pauli 
equation. Lorentz boosts from the particle rest frame point uniquely to the square-root Hamiltonian, 
but these don't exist for a massless particle. Instead, Maxwell's equations are dissected in spatial Fourier 
transform to separate nondynamical longitudinal from dynamical transverse field degrees of freedom. Upon 
their decoupling in the absence of sources, the transverse field components are seen to obey two identical 
time-dependent Schrodinger equations (owing to two linear polarizations), which have the massless free- 
particle diagonalized square-root Hamiltonian. Those fields are readily modified to conform to the attributes 
of solitary-photon wave functions. The wave functions' relations to the potentials in radiation gauge are 
also worked out. The exercise is then repeated without the considerable benefit of the spatial Fourier 
transform. 
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Introduction 



It was pointed out in a previous paper [T] that a solitary free relativistic nonzero-mass particle is described 
without any pathology whatsoever by the natural time-dependent Schrodinger equation, 

ihd{\ip{t)))/dt = ^™ 2 c 4 + |cp| 2 (i) 

whereas the widely used relativistic free-particle Klein-Gordon and Dirac equations are generally acknowledged 
not to be up to this simple task [2]. For example, negative energy solutions of the Klein-Gordon equation fail 
to be orthogonal to their positive energy counterparts that have the same momentum. This contradicts a 
fundamental property of quantum theory that makes its probability interpretation possible; unsurprisingly 
it is well-known that Klein-Gordon theory can yield negative probabilities [2]. This particular pathology of 
the second-order in time Klein-Gordon equation cannot arise if the solitary free particle is described by a 
standard first-order in time Schrodinger equation with a Hermitian Hamiltonian operator, such as that of 
Eq. (1). The particular Hamiltonian operator of Eq. (1), namely -y/m 2 c 4 + |cp| 2 , has the additional virtue of 
fully adhering to the classical Correspondence Principle, being that it is the direct quantization of the correct 
classical Hamiltonian for a solitary relativistic free particle of mass m. It is to be noted in particular that 
solitary free relativistic particles cannot have negative energies if solitary free nonrelativistic particles are to 
be restricted to having only nonnegative kinetic energies! 

In light of the problems the second-order in time Klein-Gordon equation has in describing the solitary 
relativistic free particle, Dirac appreciated the need for elementary relativistic quantum mechanics to return to 
the standard first-order in time Schrodinger equation format with Hermitian Hamiltonian operator that serves 
elementary nonrelativistic solitary-particle quantum mechanics so admirably. Lamentably, however, Dirac was 
less responsive to the exacting requirements of the classical Correspondence Principle than he was, like Klein, 
Gordon and Schrodinger, misdirectedly concerned about the fact that the solitary free particle Hamiltonian 
operator */ m 2 c 4 + |cp| 2 turns out to be a nonlocal entity in configuration representation: it doesn't seem 
to have occurred to these pioneers that this fact in no way stymies the fruitful application of perturbation 
approximations — the relativistic corrections to the atomic physics in which they were interested are obviously 
very well-suited to this approach, being compatibly small. Dirac unfortunately rejected the Correspondence 
Principle appropriate square-root Hamiltonian operator of Eq. (1) in favor of a misconceived linearization of 
it in terms of the components of the momentum operator p and the mass m, for which he argued on the basis 
of a fundamental misapprehension of the manner in which the solitary-particle time-dependent Schrodinger 
equation in configuration representation, 

ihd({v\i,{t)))/dt = (r\HW))> 

is related to the covariance requirements of special relativity [3], [4J [2] . It is clear that the operator d/dt 
on the left-hand side of this equation is the time component of the four-vector operator cd/dx^, and the 
operator H on the right-hand side of this equation is the time component of the four-vector operator cp^, 

where p^ = (H/c, p). Moreover, it was postulated by Schrodinger, and is a basic consequence of Dirac's own 
canonical commutation rule, that, 

-tW r «i#(i))) = (r|p|V(i)), 

so that the full four-vector equation, 

ihd({T\i>(t)))/dx„ = (r\Pm)), 

is guaranteed to hold in solitary-particle quantum mechanics! Since the operator d/dx^ is patently a Lorentz 
covariant four-vector, the covariance requirements of special relativity are met in solitary-particle quantum 
mechanics by simply requiring that the Hamiltonian operator H be selected so as to ensure that the four- vector 
operator p^ = (H/c, p) also transforms between inertial frames as a Lorentz covariant four- vector. This 
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requirement is automatically fulfilled by scrupulous adherence to the strongest possible form of the classical 
Correspondence Principle, i.e., that H be the quantization of precisely that classical Hamiltonian H which has 
been carefully checked to be appropriate to fully relativistic solitary-particle classical mechanics! For the free 
solitary particle of nonzero mass m, this physically methodical and highly conservative approach leaves us with 
no option but to accept Eq. (1) as its correct time-dependent Schrodinger equation description! This even 
extends to free spin \ particles of nonzero mass: notwithstanding that spin ± itself is a nonclassical attribute, 
the nonrelativistic Pauli Hamiltonian for such a particle automatically reduces to the usual nonrelativistic 
purely kinetic- energy Hamiltonian in the free-particle limit, and one can always find an inertial frame of 
reference in which a free particle of nonzero mass is completely nonrelativisticl 

Dirac, however, was much too focused on trying to cobble up a relativistic solitary-particle Hamiltonian 
operator which is local in configuration representation to be in any frame of mind to appreciate this profound 
link between the strongest form of the classical Correspondence Principle and the requirement of Lorentz 
covariance in solitary-particle quantum mechanics. Instead of pondering the details of how the requirement of 
Lorentz covariance actually impacts the time-dependent solitary-particle Schrodinger equation in configuration 
representation, Dirac was content to assume that relativistic covariance merely requires that there be essentially 
complete symmetry in the formal treatment of space and time coordinates [21 S] ■ As a result, he completely 
missed the point that the time-dependent Schrodinger equation relates the time derivative operator to an 
energy operator, neither of which are Lorentz scalars, but each of which is the time component of a Lorentz 
four-vector. Not having assimilated these basics, he conjured from whole cloth a nonexistent Lorentz scalar 
which he perceived this equation to split into two nonscalar fragments for the express purpose of displaying the 
fragment which is proportional to the time derivative on the left-hand side of the equality signl Thus primed, 
Dirac "concluded" that his phantom scalar's "completion fragment" , which is everything on the right-hand 
side of the equality sign, "must" therefore be linear in the space gradient, which suited his purpose perfectly, as 
it results in a local "Hamiltonian operator" in configuration representation! Following this "eureka moment" , 
which was the fruit of his mistakenly identifying as a scalar the time component of a four-vector, Dirac failed 
to reflect on whether a Hamiltonian operator that is linear in the space gradient, and thus in the momentum, 
could really be relativistically correct in light of the firmly established understanding that a solitary free 
particle's Hamiltonian is ineluctably the time- component of a Lorentz- covariant four-vector whose remaining 
three components are c times that free particle's three-momentum! This understanding, conjoined with the 
Lorentz transformation itself, in fact determines that the square-root Hamiltonian which occurs in Eq. (1) is 
the only correct one for the nonzero-mass free particle! 

Dirac also paid no heed to the fact that a solitary free-particle Hamiltonian operator which is linear in the 
space gradient, and therefore in the momentum operator, has, in light of Heisenberg's equation of motion, the 
unavoidable consequence that the free particle's velocity is completely independent of its momentum, which is 
an astounding contradiction of nonrelativistic free-particle physics, quantum or classical! Dirac determined the 
coefficients of his misconceived linearized Hamiltonian operator by requiring that its square be equal to the 
square of the square- root Hamiltonian operator of Eq. (1), which is a perilously weakened surrogate for the 
strong classical Correspondence Principle that produces the square-root Hamiltonian operator of Eq. (1) in the 
first place! It results in these coefficients satisfying the well-known Dirac-matrix anticommutation relations [2]. 
The free-particle velocity operator, which involves only these coefficients, is thereupon determined to equal the 
speed of light c times a three- vector of the Dirac matrices, which each square to unity. Therefore the speed of 
any free Dirac particle turns out to have the universal superluminal value V3c, irrespective of its momentum\ 
The free-particle Dirac equation in fact yields more such inordinately unphysical results. Upon using the 
misconceived linearized Dirac Hamiltonian operator in conjunction with Heisenberg's equation of motion to 
calculate the free particle's spontaneous acceleration, one finds that its magnitude has a minimum value of order 
of the "Compton acceleration" mc 3 /%, which, for the electron, is about 10 28 <7, an absolutely staggering violation 
of Newton's First Law of Motion for a free particle! The natural square-root Hamiltonian operator of Eq. (1) 
gives nil spontaneous acceleration, a result that is, of course, in complete agreement with Newton's First Law 
of Motion for a free particle. It as well gives the correct expression for the relativistic free particle's velocity in 
terms of its momentum. The extreme disparity of the results of the natural square-root Hamiltonian operator 
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versus those of the misbegotten linearized Dirac Hamiltonian operator is an overwhelming object lesson on the 
dangers inherent in any weakening of the strongest sensible version of the classical Correspondence Principle. 

Although it is routinely claimed that the Dirac equation reduces to the nonrelativistic Pauli equation for 
a spin i particle when the particle's momentum magnitude is much less than rac [2J, the "proof" of this 
assertion is definitely invalid, being unsalvageably dependent on the lapse of forgetting that at p = the two 
lower components of the Dirac spinor in the standard representation have the time-dependence phase factor 
e +(imc /h)t^ which j s totally different from the analogous factor e~^ mc that occurs in this spinor's two upper 
components! Nor does this assertion remotely accord with some of the most elementary "physics" consequences 
of the free Dirac and Pauli theories at vanishing momentum. The latter's Hamiltonian operator is just the 
nonrelativistic kinetic energy operator |p| 2 /(2m), and its particle speed operator is, of course, |p|/m. Thus a 
free Pauli particle eigenstate of vanishing momentum has vanishing speed. For the free Dirac theory, we have 
already seen that the particle speed operator is even simpler, namely the identity operator times the universal 
superluminal speed y/3 c ! So a free Dirac particle eigenstate of vanishing momentum still has this problematic 
extreme speed! 

Furthermore, notwithstanding its spin \ attribute, the free Pauli particle's orbital angular momentum is 
exactly conserved, i.e., the rate of change of its orbital angular momentum vanishes identically. The free Dirac 
particle's orbital angular momentum, however, is coupled with astonishing strength to its spin: as the free 
Dirac particle's momentum magnitude tends toward zero, the dimensionless ratio of the magnitude of the 
rate of change of its orbital angular momentum to its kinetic energy increases monotonically without bound, 
beginning from the asymptotic ultrarelativistic dimensionless ratio value v2 ! In other words, far from having 
the exactly conserved orbital angular momentum of the free Pauli particle, the free Dirac particle's spin-orbit 
torque magnitude always well exceeds that particle's own kinetic energy, and the dimensionless ratio of these 
two quantities in fact becomes arbitrarily large at low enough particle momentum! 

All of these stunningly unphysical properties of the Klein-Gordon and Dirac equations with regard to the 
description of a solitary relativistic free particle notwithstanding, and even in spite of the fact that the simple 
square-root Hamiltonian operator of Eq. (1) — which is the unique consequence of the classical Correspondence 
Principle for a solitary relativistic free particle — doesn't partake of any such pathologies, it still has always 
been the Dirac and Klein-Gordon equations, rather than Eq. (1), that are inducted into relativistic quantum 
field theory. The reason for this, of course, is that antiparticles were first observed sometime after Dirac 
began to speculate about mechanisms which could serve to hide the physically problematic unbounded-below 
negative energy spectra that are a feature of his and the Klein-Gordon equations, but which simply do not 
occur for Eq. (1). Had Dirac not been so speculating, the existence of antiparticles would have been regarded 
as an energy degeneracy of nature's full field theoretic Hamiltonian, and an explanation for that degeneracy 
would have been sought, following the grand tradition established by Wigner, Weyl and others, entirely in 
terms of the effect of a symmetry possessed by that full Hamiltonian. A particle and its antiparticle are 
distinguishable, and each can have only positive energy, so it is entirely natural that they should be described 
by two entirely independent quantum fields, with each having purely positive energy. In other words, had the 
Klein-Gordon and Dirac equations, with their problematic unbounded-below negative energy spectra never 
been concocted, it would have been perfectly straightforward to accommodate the discovery of antiparticles in 
a simple, logical framework that is very strongly grounded in physical precedent. The straightforward use of the 
purely positive energy Eq. (1) in conjunction with symmetry postulates to accommodate antiparticles has the 
theoretical advantage that it as well automatically accommodates a sensible theory of a solitary free relativistic 
particle, which the Dirac and Klein-Gordon equations are utterly unsuited to do. There is no physical reason 
whatsoever that nonrelativistic solitary particle theory should not link to relativistic particle physics in a 
completely smooth fashion, which is what Eq. (1) transparently enables. Furthermore, both the Klein-Gordon 
and Dirac equations historically arose as eccentric offshoots of Eq. (1), motivated not by legitimate physics 
concerns, but by an irrational distaste for the nonlocal character of Eq. (1) in configuration representation. 
This means that the Klein-Gordon and Dirac equations were not designed ab initio to accommodate both 
a particle and its antiparticle: this is a role into which historical happenstance has pushed them — by their 
actual patrimony they were designed to accommodate only a single type of particlel Nowadays, it is known 
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that particle-antiparticle symmetry can be slightly broken, as CP noninvariance experiments have shown 
(given the dominance of particles over antiparticles in our immediate surroundings, it would be astonishing if 
particle-antiparticle symmetry were not in fact broken). But the Dirac and Klein-Gordon fields, not having 
been designed to accommodate two particles, are highly stressed to accommodate two slightly nondegenerate 
particles, which is what corresponds to the existent symmetry breaking. It is obvious that the model with two 
independent positive-energy fields for particle and antiparticle offers vastly more flexibility to accommodate 
symmetry breaking than do the claustrophobic Dirac and Klein-Gordon models, which shoehorn two particles 
into a field structure that was designed to accommodate just one. As one example, two independent fields easily 
accommodate two slightly different masses: there is simply no way to have a single Dirac or Klein-Gordon 
field with more than one mass. 

For a nonzero-mass solitary free particle, the relativistic square-root Hamiltonian operator of Eq. (I) 
is completely determined by the Lorentz transformation. This is because one can always find an inertial 
frame in which a solitary free particle of mass m is at rest, i.e., has four-momentum (mc, 0). The Lorentz 
transformation to the inertial frame in which this particle has velocity v, where |v| < c, then takes the particle's 
four-momentum to, 

Ml - |v| 2 /c 2 )-*, mv(l - |v| 2 /c 2 r*) = (25(v)/c, p(v)), 
which, together with the identity, 

mc 2 (l - |v| 2 /c 2 )-3 = y /m?c i + |cmv| 2 (I- |v| 2 /c 2 )-\ 

implies that, 

E(\) = ^m 2 c i + |cp(v)| 2 . 

Since the classical precursor of the square-root Hamiltonian operator for the solitary free-particle of mass m 
that occurs in Eq. (1) is thus mandated by the very nature of the Lorentz transformation, it is little wonder that 
Dirac's misconceived effort to linearize the square-root character of this Hamiltonian operator has consequences 
which terribly violate well-known relativistic properties of a free particle: we have seen that these consequences 
include the blatantly unphysical universal superluminal free particle speed ^/3c irrespective of the particle 
momentum, a minimum spontaneous free-particle acceleration magnitude of order of the Compton acceleration 
mc 3 /h, namely about 10 28 <7 for the electron, which staggeringly violates Newton's First Law of Motion for a 
free particle, and the gross failure to conserve free-particle orbital angular momentum, which the nonrelativistic 
spin | Pauli theory free particle definitely does. 

For a zero-mass free particle, however, there is no inertial frame in which that particle is at rest, so we 
cannot readily derive its Hamiltonian from the Lorentz transformation, as we have done for the nonzero-mass 
free particle. Confirmation that the Hamiltonian operator given by Eq. (I) continues to be correct for a 
massless solitary free particle must be sought elsewhere. We therefore turn to the study of electromagnetic 
radiation, which is supposed to consist of massless photons. Surprisingly, we shall see that Maxwell's classical 
equations for pure electromagnetic radiation can be recast into a form that is in essence that of the time- 
dependent Schrodinger equation of Eq. (1) with m = 0. Because of the particle's vanishing mass, Planck's 
constant % can be factored out of both sides of Eq. (I), since p = —ihW r in configuration representation and 
p = hk in Fourier vector variable k-representation. That Planck's constant drops out of the relativistic solitary 
free-particle time- dependent Schrodinger equation in the m = case is a key factor in allowing that equation 
to be related to the putatively "classical" Maxwell theory. 

To make further progress, we must dissect Maxwell's four equations themselves: these are a mixed bag 
of dynamical field equations of motion and nondynamical constraint conditions. Since the time-dependent 
Schrodinger equation is purely dynamical in character, it will be necessary to properly resolve the nondynamical 
constraint conditions, a task which we now undertake. 



■5 



The electromagnetic field as a constrained dynamical system 

Since any time-dependent Schrodinger equation is linear and homogeneous, only the source-free (i.e., pure 
radiation) version of Maxwell's equations could possibly correspond to such an equation. But the resolution 
of the two nondynamical constraints amongst the the four Maxwell equations can be carried out even in the 
presence of the source terms, so we shall initially retain those terms. The four Maxwell equations for the 
electromagnetic field (E,B) with four-current source (/?,j/c) are comprised of Coulomb's law, 

V • E = p, (2a) 

Faraday's law, 

VxE = -B/c, (2b) 

Gauss' law, 

V-B = 0, (2c) 

and Maxwell's law, 

V x B = (j + E)/c, (2d) 

which, together with Coulomb's law, implies the current conservation condition, 

V-j + p = 0. (2e) 

Coulomb's and Gauss' laws both involve no time derivatives of the electromagnetic field, so they are in the 
nature of nondynamical constraints on that field, whereas Faraday's and Maxwell's law's, which both do 
involve first time derivatives of the electromagnetic field, have the character of dynamical equations of motion 
of that field. If one is presented with a set of N variables which are subject to both nondynamical equations 
of constraint and dynamical equations of motion, it is standard practice to search for N functions of those 
N variables with the property that all the dynamical equations of motion involve only a subset of N — k of 
these functions, while all the nondynamical equations of constraint involve only the remaining subset of k 
functions. The first set of N — k functions is not subject to any nondynamical equations of constraint (these 
apply exclusively to the second set of k functions) , and are regarded as a set of purely dynamical variables for 
the system. The second set of k functions, to which no dynamical equations of motion apply, may analogously 
be regarded as a set of purely nondynamical variables for the system. The equations of motion satisfied by 
the N — k unconstrained dynamical variables are then typically summarized by means of a Lagrangian or 
Hamiltonian from which they follow, respectively, via the Euler-Lagrange or classical Hamiltonian equations 
equations of motion. Hamiltonization of such a maximal set of unconstrained dynamical variables opens the 
way to the system's quantization via either the Hamiltonian phase-space path integral [5], or, equivalently, the 
slightly strengthened self-consistent extension of Dirac's canonical commutation rule [5]. 

For the electromagnetic field, N is formally infinite, but we can still usefully discuss the number of field 
degrees of freedom; e.g., the electromagnetic field (E(r, t), B(r, t)) has six field degrees of freedom. Both 
the equations of motion and those of constraint are linear for the electromagnetic field, so one can expect the 
extraction of a maximal subset of unconstrained dynamical variables (actually unconstrained dynamical fields) 
to involve appropriate linear transformations of components of the electromagnetic field (E(r, t), B(r, t)). 
Furthermore, consideration of the Coulomb and Gauss equations of constraint quickly makes it clear that 
V • E(r, t) and V ■ B(r, t) (which vanishes identically!) are each purely nondynamical single field degrees of 
freedom, and that no additional purely nondynamical field degrees of freedom are available to be extracted 
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from the six field degrees of freedom of the electromagnetic field system (E(r, t), B(r, t)). Therefore the 
electromagnetic field must have four unconstrained, purely dynamical held degrees of freedom. To cleanly 
separate the nondynamical V • E(r, t) and V • B(r, t) from the purely dynamical part of (E(r, t), B(r, t)), 
a hypothetical coordinate system in which one of the three components of the electric field E(r, t) is just 
V • E(r, t) and also in which one of the three components of the magnetic field B(r, t) is just V • B(r, t) would 
be very convenient. In such a hypothetical coordinate system, the set of the remaining two components of 
E(r, t), together with the remaining two components of B(r, t), would comprise the four unconstrained, purely 
dynamical electromagnetic field degrees of freedom. It turns out to be technically most straightforward to 
precede attempted implementation of this sort of idea by spatial Fourier transformation of the electromagnetic 
field (E(r, t),B(r, t)) and its four-current source (p(r, t), j(r, t)/c). We define, 

(E(k,t),B(M)) (27r)-i/d 3 r e - 4k - r (E(r,t),B(r,t)), (3a) 

which is the "unitary" Fourier transform. Also, 

(p(k,i)J(k,t)/c) d = f (27r)-i/d 3 r e -*- r (p(r,t),j(r,t)/ C ). (3b) 

In the penultimate section of this paper, we shall essay the trickier task of attempting to reveal the dynamical 
time-dependent Schrodinger equation character of Maxwell's source-free equations directly in configuration 
representation, without resort to this spatial Fourier transformation. It is worth remarking at this stage that 
since the charge density p(r, t) is a real-valued function, and the same is true of all the Cartesian components of 
(E(r, t), B(r, t)) and j(r, t), the corresponding spatial Fourier transforms of all these entities have the property 
that their complex conjugation is equivalent to reversing the sign of their Fourier vector argument k. Some 
key manipulations that are carried out further on rely heavily on this technical point. The spatial Fourier 
transformation of V • E(r, t) comes out be ik ■ E(k, t), which, in a coordinate system that has u^(k) = f k/|k| 
as one of its three orthogonal unit vectors, is equal to i|k| times the Uz,(k)-component of E(k, t), which we 
denote as £?i,(k, t). Coulomb's law thus obviously implies the nondynamical equation, 

E L (k,t) = -ip(k,t)/m, (4a) 

and, analogously, Gauss' law implies the nondynamical equation 

Bi(k,t)=0. (4b) 

We can therefore be quite confident that £x(k, t) and i?L(k, t) exhaust the nondynamical components of E(k, t) 
and B(k, i) respectively, and that the remaining two components of each of these two fields will be purely 
dynamical, i.e., free of any nondynamical constraint. But to demonstrate this in detail, we must explicitly 
display the remaining two mutually orthogonal unit vectors, which are each as well orthogonal to Uz,(k), and 
then work out the consequences of the Maxwell equations for the remaining two components of both E(k, t) 
and B(k, £)) in that coordinate system, in order to verify that purely dynamical equations of motion which 
involve only these four components result. In the particular case that the four-current source (p(k, t), j(k, t)/c) 
vanishes, we also need to demonstrate that the now homogeneous equations of motion obtained for these four 
unconstrained dynamical components of (E(k, t), B(k, t)) are equivalent to the schematic Schrodinger Eq. (1) 
with m = — note as well that in this source-free case the two nondynamical components El(\£, t) and i?i(k, t) 
of (E(k, t), B(k, t)) vanish identically, as is seen from Eqs. (4). 

In order to obtain two mutually orthogonal unit vectors which are both also orthogonal to the unit vector 
Ui(k) = k/|k|, we display u^(k) in Cartesian coordinates: it is simply the well-known unit vector in the radial 
direction that the Fourier vector argument k points toward, expressed in terms of that vector's spherical polar 
angles </>k and #k, 
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u L (k) = (cos 0k sin k , sin k sin k , cos 0k). (5a) 

Now because u £ (k) = k/|k|, u £ (— k) = — u L (k), i.e., u^(k) has the same odd parity that k has. Therefore 
the parity flip mapping k — > k corresponds to the polar angular mapping 0k — > 0k + tt, because this sends 
sin 0k — ► — sin 0k and cos 0k — ► — cos 0k, thus sending, from Eq. (5a), u^(k) — »■ — u^(k). However, if we instead 
choose to carry out the polar angular mapping 0k — > 0k + tt/2, then sin 0k — > cos 0k, cos 0k — > — sin 0k, and 
Ui(k) — > Ui(k), where, 

ui(k) = (cos </>k cos 0k, sin 0k cos 0k, — sin 0k), (5b) 

is readily checked to be a unit vector that is orthogonal to Ui(k). The parity flip angular mapping 0k — > 0k + 7r 
reveals that Ui(k) is ako of oc?c? parity. With the mutually orthogonal unit vectors u^(k) and ui(k) in hand, 
we can now readily construct a third unit vector u 2 (k) which is orthogonal to both of these, 

u 2 (k) = u L (k) x ui(k) = (-sin k , cos k , 0). (5c) 

It is immediately seen that u 2 (k) is of even parity. By using the identity a x (b x c) = b(a • c) — c(a • b), 
or, alternatively, the spherical polar angular representations given by Eqs. (5a)-(5c), it is readily checked that 
Ui(k), u 2 (k) and u^(k) comprise a "right-handed" orthonormal local vector triad, i.e., 

ui(k) x u 2 (k) = u L (k), u 2 (k) x u L (k) = m(k) and u L (k) x m(k) = u 2 (k). (5d) 

Turning now to the implications of Maxwell's equations in this coordinate system, we have already noted 
that the Coulomb and Gauss laws imply the two nondynamical Eqs. (4a) and (4b). Upon spatial Fourier 
transformation, Faraday's law, Eq. (2b), becomes, 

ik x E(k, t) = -B(k, t) jc. (6a) 

Noting that k = |k|ui,(k), and that, 

E(k, t) = E x (k, t)ui (k) + E 2 (k, t)u 2 (k) + E L (k, t)u L (k) , 

where, 

£i(k,f) d = m(k) • E(k,t), E 2 (k,t) d = u 2 (k) • E(k,f) and E L (k,t) d = u L (k) • E(k,i), 

and analogously for B(k, t), for which the Gauss law result embodied by Eq. (4b) already permits us to conclude 
that £?£,(k, t) = 0, we apply Eq. (5d) to the left-hand side of Eq. (6a), and thereby obtain the two additional 
equations, 

*Bi(k,t) = -|ck|fb(k,t), (6b) 

and, 

iB 2 (M) = |ck|f?i(k,t). (6c) 

Before we turn to Maxwell's law, Eq. (2d), it is convenient to treat the current conservation condition, Eq. (2e), 
which is a constraint on the four-current source that follows from Maxwell's and Coulomb's laws. Upon spatial 
Fourier transformation, Eq. (2e) becomes, 

ik-j(k,t) = -p(k,i), (6d) 
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which immediately yields the longitudinal source current component in terms of the rate of change of the 
charge density, 

j L {k,t)=ip(k,t)/\k\. (6e) 

Upon spatial Fourier transformation, Maxwell's law, Eq. (2d), becomes, 

ik x B(k, t) = (j(k, t) + E(k, t))/c. (6f) 

The left-hand side of Eq. (6f) has a vanishing component in the Ui(k)-direction, and the joint consequence of 
that and the Coulomb law result embodied by Eq. (4a) for its right-hand side is simply the constraint on the 
longitudinal source current component that is embodied by eq. (6e). More interesting are the two equations 
that follow from the components of Eq. (6f) in the u x (k) and u 2 (k) directions — these bear a strong resemblance 
to the Eqs. (6b) and (6c) which follow from Faraday's law, 

i£?i(k,t) = \ck\B 2 {k,t) -iji(k,t), (6g) 

and, 

iE 2 {k,t) = -|ck|Bi(k,t) -y 2 (k,i). (6h) 

Aside from the purely source constraint requirement of Eq. (6e) and the reconfirmation that Bl (k, t) must 
vanish, which is already a consequence of Eq. (4b) (which is itself the result of the Gauss law), the Faraday 
and Maxwell laws have yielded four dynamical equations of motion, namely Eqs. (6b), (6c), (6g) and (6h), 
which involve only the four transverse field components E\{k, t), B 2 (k, t), E 2 (k, t) and £?i(k, t). Absolutely no 
nondynamical equations of constraint for any of these four transverse field components have eventuated from 
any of the Maxwell equations. It is therefore clear that the six field degrees of freedom of (E(k, t), B(k, t)) 
have now been successfully partitioned into four unconstrained, purely dynamical transverse field degrees of 
freedom and two purely nondynamical longitudinal field degrees of freedom E L (k, t) and B L (k, t), whose values 
are actually given by the simple nondynamical constraints of Eqs. (4a) and (4b). In addition, it has, of course, 
transpired that the four-current source (p(k, t), j(k, t)/c) cannot be chosen arbitrarily, but is subject to the 
source constraint given by Eq. (6e) . 

Linear algebraic decoupling of the four transverse dynamical fields 

Eq. (6g) for the dynamical transverse fields E\(k,t) and B 2 (k, t) is clearly coupled to Eq. (6c), and likewise 
Eq. (6h) for the dynamical transverse fields E 2 (k,t) and Bi(k,t) is clearly coupled to Eq. (6b). Some investi- 
gators may be tempted to decouple these equations by taking second time derivatives, but such an approach is 
entirely unnecessary and involves a risk of introducing extraneous solutions that don't actually apply to these 
equations — indeed taking an unwarranted second time derivative is precisely how the unphysical, unbounded- 
below negative energy spectrum was inadvertently forced into the relativistic Klein-Gordon equation for a 
nonzero-mass free particle. Eqs. (6g) and (6c) are easily decoupled by the straightforward expedient of taking 
their sum and difference, and the same applies to Eqs. (6h) and (6b). Adding Eq. (6c) to Eq. (6g) yields, 

idiE^t) + B 2 (k,t))/dt = |ck|(Ei(k,i) + B 2 {k,t)) - iji(k,t), (7a) 

while subtracting Eq. (6b) from Eq. (6h) yields, 

id(E 2 (k,t) - B 1 (k,t))/dt = \ck\(E 2 (k,t) - Si(k,t)) - y 2 (M). (7b) 
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One can also subtract Eq. (6c) from Eq. (6g) to obtain, 

i0(£i(k,i) - B 2 (k,t))/dt = -IckK^M) - B 2 (M)) - yi(k,t), (7c) 
and add Eq. (6b) to Eq. (6h) to obtain, 

id(E 2 (k,t) + B 1 (k,t))/dt = -\ck\(E 2 (k,t) + B 1 (k,t)) - y 2 (k,t). (7d) 

Now it turns out that Eq. (7c) is not independent of Eq. (7a); in fact, Eq. (7c) is actually equivalent to 
Eq. (7a)! The reason for this is somewhat involved; it is related to the previously mentioned fact that for any 
Cartesian component of E(k,t), B(k,t), or j(k,t), complex conjugation is equivalent to changing the sign of 
the Fourier vector argument k. Making matters a bit more complicated is the fact that E\ (k, t) and E 2 (k, t) 
are not Cartesian components of E(k,t) because £?i(k, t) = E(k,i) • Ui(k) and E 2 (k,t) = E(k,t) • u 2 (k). 
Because ui(k) is of odd parity in its argument k, complex conjugation of Ei(k,t) is equivalent to changing 
both the sign of its argument k and its overall sign! However, because u 2 (k) is of even parity in its argument 
k, complex conjugation of E 2 (k,t) is equivalent to merely changing the sign of its argument k. Exactly the 
same distinction with regard to complex conjugation holds between B\(k, t) and B 2 (k, t), and as well between 
ji(k, t) and j 2 (k, t). Now if we take the complex conjugate of both sides of Eq. (7c) and apply what we have 
just learned, the result is, 

-*0(-£i(-k,t) - B 2 (-k,t))/8t = -|ck|(-f?i(-k,t) - B 2 (-M)) - iji(-M). (7c) 

Upon combining signs in Eq. (7e), we find that it resembles Eq. (7a) in every respect, except for the fact that 
all occurrences of the Fourier vector argument k have effectively had their sign reversed. However, because 
Eq. (7c) is supposed to hold irrespective of what value is assumed by k, we are free to make the simple one-to- 
one formal transformation k — > k which turns Eq. (7e) into Eq. (7a). Furthermore, if we take the complex 
conjugate of both sides of Eq. (7d) and apply to it what we have learned above, the result is, 

-id(E 2 (-k,t) - B!(-k,t))/0i = -\ck\(E 2 (-k,t) - Bi(-k,i)) + y 2 (-k,t). (7f) 

Upon reversing the sign of both sides of Eq. (7f), we find that it resembles Eq. (7b) in every respect, except 
for the fact that all occurrences of the Fourier vector argument k have effectively had their sign reversed. But 
we are, of course, again justified in making the simple one-to-one formal transformation k — > k which turns 
Eq. (7f) into Eq. (7b). 

We have thus succeeded in replacing the four coupled equations of motion for the dynamical transverse 
fields by two nontrivially complex- valued and fully decoupled such equations, namely Eqs. (7a) and (7b). If we 
multiply both of these equations through by K, and then set both of the transverse source currents ji(k, t) and 
j 2 (k,t) to zero, Eqs. (7a) and (7b) assume precisely the schematic form of Eq. (1) with m — 0, i.e., they are 
of the form of time-dependent Schrddinger equations for a solitary relativistic massless free particle. The fact 
that there are two such equations suggests, in light of the detailed electromagnetic field composition of each 
of their two apparent "wave functions" , that they describe the amplitudes for two linear polarization states of 
the solitary massless particle. We shall now further investigate this interesting source-free limit of Maxwell's 
equations. 

The Schrodinger character of the source-free Maxwell equations 

When the four-current source (p(k, t), j(k, t)/c) vanishes altogether, Eqs. (4) show that the two nondynamical 
longitudinal electromagnetic field elements £x(k, i) and Bj,(k,i) vanish identically as well. The only physics 
that remains is purely dynamical and transverse, and is fully dcscriblcd by the two relativistic, massless, solitary 



10 



free-particle Schrddinger-style equations, 

ihd(E!(k,t) +B 2 (k,t))/dt = \chk\(E 1 (k,t)+B 2 (k,t)), (8a) 

and, 

ihd(E 2 (k,t) - Bi(k,t))/0t = \chk\(E 2 (k,t) - Bi(k,t)), (8b) 

which follow from Eqs. (7a) and (7b) in the source-free situation. The detailed structure of the two putative 
"wave functions" in terms of the transverse electromagnetic field components strongly suggests that they 
represent the amplitudes for the two possible transverse linear polarization states of the solitary, massless, 
free electromagnetic field particle. There is a technical snag, however, which bars such an interpretation from 
being immediately made: the "wave functions" that appear in the Schrddinger equations of Eqs. (8) are sums 
and differences of transverse electromagnetic field components, which have the character of energy- density 
amplitudes, whereas true solitary-particle wave functions have the character of probability- density amplitudes. 
To get a feeling for just what energy is represented by the two complex- valued "wave functions" of Eqs. (8), we 
wish to integrate the sum of their absolute squares over all of the Fourier vector- variable k-space. We begin 
by integrating over just the absolute square of the "wave function" of Eq. (8a), 

ftPkQEh. + B 2 \ 2 = Jd 3 k(\E 1 \ 2 + \B 2 \ 2 ) + ja?k{E{B 2 + B 2 E{), (9a) 

where we have temporarily suppressed writing out the arguments of the transverse field components to save 
space. However, bearing in mind the discussion between Eqs. (7d) and (7e), we have that, (E\{k,t))* — 
-.Ei(-k,i) and B 2 {k,t) = (B 2 (-k,t))* , from which we readily deduce that f d 3 kE$B 2 = -/d 3 kB 2 *£i, and 
therefore that the second integral on the right-hand side of Eq. (9a) vanishes. Analogous arguments show 
that when one integrates over the absolute square of the "wave function" of Eq. (8b), the integration over 
the corresponding two cross terms vanishes as well. Therefore, the result of integrating over the sum of the 
absolute squares of these two "wave functions" is, 

/^kflEi + B 2 \ 2 + \E 2 - Bil 2 ) = fd?k(\E 1 \ 2 + \E 2 \ 2 + \Bi\ 2 + \B 2 \ 2 ). (9b) 

Now the integral on the right-hand side of Eq. (9b) is equal to twice the total energy present in the transverse 
electromagnetic field components, which is, of course twice the total electromagnetic field energy, since the 
longitudinal components of the electromagnetic field vanish identically in the source-free case. Therefore the 
energy represented by the two complex- valued "wave functions" of Eqs. (8) is twice the total electromagnetic 
field energy. 

Now let us suppose that the sole source of the electromagnetic field energy which is present is a solitary 
photon. That photon's energy density in k-space (which is effectively momentum-space, since tik is the pho- 
ton's momentum) is then equal to half of the sum of the absolute squares of the two complex-valued "wave 
functions" of Eq. (8), as we have learned from Eq. (9b). Now insofar as the solitary photon has its Fourier 
vector variable equal to k, i.e., insofar as it has momentum Tik, it clearly also has energy |c?ik|. Therefore, 
we can convert our photon's energy density in k-space — which is half of the sum of the squares of the two 
complex-valued "wave functions" of Eqs. (8) — to its probability density in k-space by simply dividing that 
energy density by |c?ik|. This implies that we can convert each of the two transverse electromagnetic field 
"wave functions" of Eqs. (8) to a proper solitary photon wave function (whose absolute square yields a proba- 
bility density) by dividing it by (2|c7ik|) 2 . It is clear that both of these proper solitary photon wave functions 
will satisfy the very same Schrddinger equation that the two transverse electromagnetic field component "wave 
functions" of Eqs. (8) satisfy: the factor of (2|c7ik|)~ 2 doesn't interfere with the validity of that time-dependent 
Schrddinger equation. Therefore, when only a solitary photon is present, its two linear polarization wave func- 
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tion components (complex- valued probability amplitudes) are given in terms of the corresponding transverse 
electromagnetic field components by, 

(k|^(i)) = (2\cfik\)-i(E 1 (k,t) + B 2 (k,t)), (10a) 

and, 

(k|^,(t)> - (2\chk\)-i(E 2 (k,t) - Bi(k,t)). (10b) 

It is convenient to as well explicitly write down the parity-reversed complex conjugates of these solitary-photon 
linear polarization wave function components, 

{Mt)\ -k) = -(2|c&k|)-i(£i(k,t) - MKt)), (10c) 

and, 

(lM*)|-k) = (2|c?ik|)-i(£ 2 (k,t) + S 1 (k,t)), (lOd) 

because, with these in hand, the relationships of the solitary-photon linear polarization wave function compo- 
nents to the transverse electromagnetic field components can be inverted, 

E x (k,t) = (|c&k|/2)+5«k|Vi(t)) - (Mt)\ -k)), (10e) 
£ 2 (M) = (|cRk|/2)+i((k|^(t)> + (Mt)\ - k)), (10f) 



iMM) = -(|cftk|/2)+*«k|^(t)) - (^(t)| - k)), (lOg) 



B 2 (k,t) = (|cfck|/2)+*«k|^i(i)) + (^i(i)l - k)). (lOh) 

It is worth explicitly reiterating that the two complex- valued linear polarization components of the solitary- 
photon wave function satisfy the massless case of the relativistic free-particle time-dependent Schrodinger 
equation that is given by Eq. (1), 

ihd((k\Mt)))/dt = \chk\ (k|Vi(t)), (lla) 



ihd{(k\Mt)))/9t = \chk\ {k\i> 2 {t)}. (lib) 

Finally, it is worthwhile to relate the solitary free-photon's complex- valued wave function to the components 
of the electromagnetic four-vector potential to which it corresponds. The electromagnetic four-vector potential 
does have a gauge ambiguity issue which unfortunately is not fully resolved by the relativistically invariant 
Lorentz condition — suppression of the ensuing timelike and longitudinal "ghost radiation" [7j [8] requires a 
further stipulation: probably the most intuitively appealing is to require the scalar potential to be uniquely 
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determined, in strictly homogeneous and causal fashion, by the charge density, which is, after all, its notional 
source after imposition of the Lorentz condition. This produces results that are no less definite than those of 
the Coulomb gauge — in fact these two gauges produce identical results for all static charge densities — but with- 
out the Coulomb gauge's disconcerting instantaneous scalar potential response at arbitrarily large distances 
to charge density change. In the present source-free case, both gauges are, in fact, identical to the radiation 
gauge 9 , V • A(r, t) — (f>(r, t) — 0, which causes the four- vector potential to have only four transverse dynam- 
ical field degrees of freedom, in complete agreement with the situation discussed above for the electromagnetic 
field in this source-free case. The relation of the electromagnetic field to the four- vector potential is, of course, 
given by, 

B(r,i) = V x A(r,i), (12a) 

and, 

E(r,t) = -V#r,t)-A(r,t)/c (12b) 
which, in spatial Fourier transform become, 

B(k,t) = ik x A(k,i), (12c) 

and, 

E(k, t) = -ik(j)(k, t) - A(k, t) /c. (12d) 

Upon applying to it Eq. (5d), Eq. (12c) readily yields the two transverse components of A(k, t) in terms of 
those of B(k, t), 

A 1 (k,t) = -iB 2 (k,t)/\k\, (13a) 



A 2 (k,t)=iB 1 (k,t)/\k\, (13b) 
and Eq. (12d) immediately yields the two transverse components of A(k, t) in terms of those of E(k, t), 

Ai(M) = -cEi(k,t), (13c) 

A 2 {k,t) = -cE 2 (k,t), (13d) 
Upon putting Eqs. (lOe) through (lOh) into Eqs. (13a) through (13d) above, we obtain, 

A 1 (k,t) = -i(|cnk|/2)+i((k^i(i)) + (Mt)\ - k»/|k|, (14a) 

A 2 {k,t) = -»(|c?ik|/2)+5«k|^(t)> - (Mt)\ - k»/|k|, (14b) 
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= -c(|cftk|/2) + *(<k|Vi(t)) - (V>i(*)l - k)), (14c) 

A 2 (k,t) = -c(|cRk|/2) + i«k|^(t)) + (Mt)\ - k ))- ( 14d ) 
Eqs. (14a) through (14d) can now be inverted, 

(MM*)) = (2|cftk|)-i(*|k|^ 1 (k,t) - ii(k,t)/c), (15a) 

(k|^ 2 (*)) = (2|c^k|)-i (i|k|A 2 (k, i) - i 2 (k, t)/c), (15b) 

(V'lWI-k) =(2|cftk|)-i(i|k|Ai(k,t)+ii(k,t)/c), (15c) 

(^(t)|-k) = -(2|cftk|)-3(i|k|A 2 (k,t) + i4 2 (k,t)/c). (15d) 



From Eqs. (15a) and (15b) it is apparent that the correct Schrodinger equation quantization of the solitary 
free photon requires not only the two transverse components of A(k, t), but as well the two transverse com- 
ponents of A(k, i). The two linear polarization state wave function components, (k|^i(i)) and (k|^> 2 (£)), are 
each ineluctably complex-valued objects in a way that is thoroughly nonsuperficiai. it requires two "classi- 
cal" field degrees of freedom, such as both Ei(k,t) and S 2 (k, t), or both A\(\s.,t) and A\(\s.,t), to comprise 
one such deeply complex- valued quantum wave function component. Of course this bodes well for the next 
level of quantization, wherein our Schrodinger equation wave function components are themselves promoted 
to become operators which have prescribed commutation relations with their own Hermitian conjugates: this 
reflects their complex-valued makeup from independent fields which are interpreted as being mutually canon- 
ically conjugate, a status for which the pair Ai(k, t) and Ai(k, t) are, of course, prime candidates. We see 
that the automatic solitary photon "first quantization" that is simply part and parcel of the very nature of 
Maxwell's supposedly "classical" equations also automatically has properties which anticipate and facilitate 
"second quantization" . Once transverse source currents are present, the solitary photon Schrodinger equation 
becomes inhomogeneous, i.e., it no longer is a Schrodinger equation, as we clearly see from Eqs. (7a) and (7b). 
The inhomogeneity of what, in the source-free case, had been the solitary photon Schrodinger equation, of 
course bespeaks the creation and destruction of such photons. It is quite remarkable, however, just how well- 
organized the solitary photon wave function is ab initio for rising to the challenges of the eventually necessary 
"second quantization" . 

Configuration space approach to the Schrodinger character of Maxwell's source- 
free equations 

Having gained insight from dissection of the Maxwell equations in spatially Fourier-transformed formulation, 
we now try our hand at teasing out the Schrodinger character of their transverse dynamical segment directly 
in configuration representation. We have learned that "dynamical" and "transverse" are effectively synonyms 
for Maxwell's equations, so we simply focus on separating the electromagnetic fields E(r, t) and B(r, t), and 
also the current j(r, t), into their physically natural transverse and longitudinal parts. This, of course, requires 
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no action whatsoever for B(r, t), which Gauss' law, Eq. (2c), marks a purely transverse. The longitudinal part 
of E(r, t) is the negative of the gradient of the very same scalar potential <fi(r,t) which describes E(r, t) fully 
in the electrostatic limit, i.e., 

Ei(r,t) = -V0(r,i), (16a) 

where, of course, 

0(r,t) (47r)-7rf 3 rV(r',i)/|r-r'|. (16b) 

Because of the Green's function identity, 

V2(l/|r-r'|) = -4^,5( 3 ) (r-r'), 

we have that, 

V 2 ^(r,t) = -p(r,t), (16c) 

and, 

V-E(r,i) = V-E L (r,f)=p(r,i), (16d) 

as required by the Coulomb law, Eq. (2a). Eq. (16d) follows from Eq. (16a), Eq. (16c) and, 

E(r,t) = E L (r,t)+Er(r,t), (16e) 

the separation of E(r, t) into its longitudinal and transverse parts, where, of course, by definition, 

V • E T (r,t) = and V x E L (r,i) = 0. (16f) 

Now the current conservation condition given by Eq. (2e) has a formal structure that is very similar to that 
of Coulomb's law, and therefore permits an analogous determination of jz,(r, t), the longitudinal part of the 
current j(r, t), 

jL(r,t)=V0(r,t), (16g) 

which, together with Eq. (16c) and the relations for j(r, t), jz,( r >*) and jr(r, t) which are analogous to Eqs. (16c) 
and (16f), implies that, 

V • j(r,t) - V ■ ji(r,t) = -p(r,t), (16h) 

as required by the current conservation condition of Eq. (2e). 

We are now in a position to reexpress the Maxwell and Faraday laws in terms of only the transverse fields 
Er(r,i) and B(r, t), and the transverse current jr(r,t). In particular, Eqs. (16g) and (16a) permit us to 
deduce that, 

j(r,t) + E(r,t) = j T (r,i) + E T (r,i), 
which permits the Maxwell law of Eq. (2d) to be rewritten, 

V x B = (j T + E T )/c. (17a) 
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From Eqs. (16a) and (16e), or from Eq. (16f), we can deduce that V x E(r, t) = V x Ey(r, t), which permits 
the Faraday law of Eq. (2b) to be rewritten, 

V x E T = -B/c. (17b) 

This completes the divorce of the dynamical Maxwell and Faraday laws from the purely nondynamical variable 
E^(r, t), whose value is given in detail by Eqs. (16a) and (16b). We can now combine the purely dynamical 
and transverse versions of the Maxwell and Faraday laws which are given by Eqs. (17a) and (17b) into a single 
complex- valued equation by adding Eq. (17b) to Eq. (17a) multiplied through by the imaginary unit i. This 
is readily seen to produce the result, 

id(E T (r,t) +iB(r,t))/dt = cVx (E T (r,t) + iB(r,t)) - ij T (r,t). (17c) 

If we instead subtract Eq. (17b) from Eq. (17a) multiplied through by the imaginary unit i, we obtain, 

id(E T (r,t) - iB(r,t))/dt = -cV x (E T (r,i) - iB(r,t)) - ij T (r,t), (17d) 

The difference between Eqs. (17c) and (17d) above bears a strong similarity to the difference between Eqs. (7a) 
and (7c) or that between Eqs. (7b) and (7d). As in those instances, it is readily shown that Eqs. (17c) and 
(17d) above are, in fact, equivalent. We kept Eqs. (7a) and (7b) in preference to Eqs. (7c) and (7d) because, 
when the transverse current source terms were dropped, the former two evidenced a manifestly nonnegative 
Hamiltonian operator, which is physically appropriate for a solitary free particle. Here, unfortunately, neither 
of the candidate Hamiltonian operators (cWx) and (—cKVx) turns out to be nonnegative, and both have the 
further peculiarity of being odd parity operators, which is unacceptable for a solitary free photon Hamiltonian 
operator. We shall, in fact, need to meld a piece from each of these two Hamiltonian operators to one another 
in such a way as produces in configuration representation the physically appropriate Hamiltonian operator for 
the solitary free photon which Eqs. (7a) and (7b), sans their transverse current source terms, have already 
delivered to us in Fourier vector variable representation as \cKk\. But before we discuss this melding of pieces 
from the two Hamiltonian operators, we need to look at the details of the passage to the source-free situation, 
and we also need to pass from the complex-valued transverse- vector electromagnetic field strengths that are 
found in Eqs. (17c) and (17d) to complex-valued transverse-vector probability density amplitudes, since it is 
these which befit the quantum description of a solitary particle (the transverse- vector character of these proba- 
bility density amplitudes reflects the the fact that the particle has only two polarization degrees of freedom). In 
the source-free case that j(r, t) = and p(r,t) — 0, it follows that (fi(r,t) — 0, and therefore also Ei,(r,i) = 0, 
which, in turn, causes E-r(r,i) to be equal to E(r, t). Thus, after they are multiplied through by h, Eqs. (17c) 
and (17d) become the homogeneous time-dependent Schrodinger-like equations, 

ihd(E(r,t) +iB(r,t))/dt = chV x (E(r, t) + iB(r, t)), (17e) 

and, 

ihd(E(r,t) - iB(r,t))/dt = -cftV x (E(r,t) - iB(r,t)), (17f) 

whose respective "wave functions" (E(r, t) ± «B(r,t)) are complex-valued vectors that are purely transverse, 
since, in the source-free case, V • E = from the Coulomb law, as well as V • B = from the Gauss law. 
Their respective "Hamiltonian operators" (±c?iVx) are Hermitian on such a "wave function" space, since it 
is readily shown, using integration by parts that, 

Jd 3 r [(a(r, t))* ■ (±V x b(r, t))} = fd*r [(±V x a(r, t))* ■ b(r, t)} . (17g) 
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Because the total purely electromagnetic field energy, which is \ J d 3 r (|E(r, t)\ 2 + |B(r, t)\ 2 ), is readily 
shown to be independent of time in the source-free case (indeed this follows from the real and imaginary parts 
of either of Eqs. (17c) or (17f), in conjunction with Eq. (17g)), we can define two candidate solitary-particle 
complex- valued vector wave functions which are properly normalized to unity, namely 

*(r,t) d = f (E(r,t) +iB(r,i))/(/dV (|E(r',t)| 2 + |B(r', t)\ 2 ))i , (18a) 

and also (\I>(r, t))*, the complex conjugate of this *ff(r,t). It is clear from Eq. (18a) that, 

Jd 3 r [(¥(r,t))*-¥(r,t)] = 1, (18b) 

a property which (H>(r,i))* obviously shares, 

/d 3 r[*(r,t)-(*(r,t))*] = l. (18c) 

We further note that, because in the source-free case V • E(r,i) = and V • B(r,t) = 0, both *(r,i) and 
(\P(r,t))* are strictly transverse, i.e., 

V • *(r, t) = and V • (*(r, t))* = 0, (18d) 

but that, as we see from Eqs. (17c) and (17f), the time-dependent Schrodinger equations which *(r,t) and 
(\E r (r, t))* obey have Hamiltonian operators of opposite sign, 

ih*(r,t) = chV x *(r,i), (18c) 



ift(*(r,t))* = -cW x (¥(r,t))*. (18f) 

Unfortunately, neither of the Hamiltonian operators (cSVx) and (-cWx) in the two Schrodinger equa- 
tions just given is nonnegative, which, in light of the nonnegativity of the purely electromagnetic field energy 
| J d 3 r (|E(r,i)| 2 + |B(r, i)| 2 ), is an absolutely necessary requirement for the physically proper description of a 
solitary free photon. In addition, both of these Hamiltonian operators are odd parity operators, which implies 
that reversing the parity of any of their energy eigenstates produces an energy eigenstate whose eigenenergy has 
the opposite sign to that of the eigenenergy of the original energy eigenstate! This complete energy eigenspec- 
trum sign symmetry is again incompatible with the need for solitary free photons to be of strictly nonnegative 
energy. The nonconservation of parity for a free solitary photon which such a Hamiltonian implies is as well 
incompatible with the conservation of parity on the part of electromagnetic theory. 

In order to overcome these problems, we must meld the nonnegative part of the eigenenergy/eigenstate 
spectrum of (cfiVx) with the complementary positive part of the eigenenergy/eigenstate spectrum of (-cWx). 
The resulting nonnegative Hamiltonian operator will clearly be the absolute value of the operator (c/iVx), i.e., 
the operator \cKS7 x |. To gain further insight into the precise nature of this melded Hamiltonian operator 
|c?iV x |, let us look at the purely transverse eigenstates of (cWx). Each of these eigenstates is clearly a 
complex- valued transverse vector field which consists of a single Fourier component whose Fourier vector vari- 
able direction is strictly orthogonal to the longitudinal unit vector m,(k). Upon operating on such a single 
Fourier component proportional to e lkr , the operator (Vx) clearly becomes i(kx) = i|k|(ui,(k)x), whereas 
the two-dimensional transverse Fourier vector variable space is clearly fully spanned by the two transverse unit 
vectors Ui(k) and U2(k), or, much more usefully in this particular instance, the two complex-valued transverse 
unit vectors (2) _ '(ui(k) ± iu 2 (k)). By making use of the cross product identities of Eq. (5d) in conjunction 
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with the facts just discussed, one readily verifies that the purely transverse eigenstates of the Hcrmitian oper- 
ator (Vx) are of the form, 

< ±} (r) ( Ul (k) ±*u 2 (k)) 3 ( ± )(k)e lk - r , (19) 

and have the corresponding eigenvalues ±|k| for the operator (Vx). Therefore the corresponding two eigenval- 
ues of the melded Hamiltonian operator \cKV x | on the above two purely transverse eigenstates are both simply 
|cfik|, in complete agreement with those of the fully diagonalized Hamiltonian for the two time-dependent 
Schrodinger equations that are given by Eqs. (11a) and (lib). Thus we clearly see that an equivalent way to 
express the action of our melded configuration-space Hamiltonian operator |c?iV x | on the purely transverse 
subspace is as simply \ch~k\, which can be loosely styled as |c(— ihV)\ in configuration representation. Because 
of the purely transverse, two-dimensional nature of the subspace on which it is obliged to operate, the non- 
negative melded configuration-space Hamiltonian |c7iV x | loses all trace of its curl character and becomes 
indistinguishable from the loosely styled \c(— iKV)\, thus falling into line with the m = case of the Eq. (1) 
Hamiltonian, namely |cp| — |c?ik| = |c(— iKV)\. 

Summarizing, we have that the solitary free-photon wave function in configuration representation is a 
complex- valued vector field *(r, t) which is strictly transverse, i.e., 

V-*(r,i)=0, (20a) 

generally normalizable to unity (unless idealized), 

Jd 3 r [(*(r,t))*-*(r,t)] = 1, (20b) 
and satisfies the time-dependent Schrodinger equation, 

ih&(r,t) = \c(-ihV)\V(r,t), (20c) 

in accord with the m = instance of Eq. (1). The loosely styled configuration-space Hamiltonian operator 
|c(— i?iV)| of the above Schrodinger equation is actually, in configuration representation, a symmetric nonlocal 
integral operator (something which Klein, Gordon, Schrodinger and Dirac sought desperately to avoid, thereby 
playing havoc with physical cogency), whose kernel is given by, 

(r| \c(-ihV)\ |r') = (r| |cftk| |r') = -(ch/(2ir 2 R))d 2 (R/(R 2 + e 2 ))/dR 2 , (21) 

def 

where R = |r — r'| and e is an infinitesimal length. 
Conclusion 

Finally, it is to be noted that the very first quantum theorist was James Clerk Maxwell. His celebrated 
equations faithfully encompassed the correct quantum description of the solitary free photon long before 
Erwin Schrodinger was to accomplish the same feat for the solitary nonrelativistic free particle. And Maxwell's 
formidable theoretical physics machinery already yielded up the first instance of intrinsic particle degrees of 
freedom, with all their subtlety! By the grace of almost cosmic coincidence, Maxwell, unlike any of his quantum 
theory successors, could accomplish all this with no use whatsoever of Planck's world-changing constant, which, 
still undiscovered, silently awaited the future — the massless nature of the photon permits Maxwell's magnificent 
equations to simply slip away from Ti's grasp. 
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